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ABSTRACT .  A  simple  formula  is  described  that  can 
extrapolate  the  kill  probability  data  of  NAVORD  Report  7019, 

(NOTS  TP  238Z),  "A  Handbook  on  the  Effectiveness  of  Cluster 
Weapons  Against  Unitary  Targets",  by  Eldon  L.  Dunn,  with  S' 
reasonable  degree  of  accuracy.  There  are  limitations  in  usiife 
the  formula,  however;  and  since  the  formula  has  an  empirical 
origin,  its  limitations  are  discussed  at  length  with  the  aid  of  ||  «  ^ 
some  theory,  several  graphs,  nomf»  rmnnprir'al  AVaTnnljiRSs  fw  pn 


and  some  numerical  examples p 
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FOREWORD 


This  report  is  intended  as  a  supplement  to  NAVORD  Report  7019, 
NOTS  xi?  2382,  entitled  "A  Handbook  on  the  Effectiveness  of  Cluster 
Weapons  Against  Unitary  Targets",  by  Eldon  L.  Dunn.  A  technique 
is  presented  which  can  be  used  to  extrapolate  kill  probability  data 
from  the  Handbook  to  determine  the  effectiveness  of  larger  numbers 
of  cluster-delivered  bomblets  and  larger  delivery  patterns. 

The  work  was  carried  out  at  the  U.  S.  Naval  Ordnance  Test  Station 
under  Bureau  of  Naval  Weapons  WepTask  RMM042-004/216-1/F008- 
08-006. 

This  report  is  released  at  the  working  level  for  informational 
purposes  only. 
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INTRODUCTION 


Recent  studies  to  determine  the  effects  of  very  large  numbers  of 
cluster-delivered  bomblets  and  the  effects  of  the  relatively  large  actual 
delivery  patterns  have  led  to  a  search  for  a  method  of  extrapolating 
data  which  already  exists  in  NAVORD  Report  7019,  "A  Handbook  on 
the  Effectiveness  of  Cluster  Weapons  Against  Unitary  Targets".^  The 
graphs  in  the  Handbook  cover  the  delivery  of  up  to  400  bomblets  in  a 
cluster  and  delivery  standard  deviations  of  up  to  8  target  widths* 


In  studying  the  "Handbook",  there  seemed  to  be  a  remarkable 
similarity  between  many  of  the  graphs.  Further  encouragement  was 
received  by  superimposing  one  graph  upon  another.  In  some  cases, 
the  graphs  were  virtually  identical  except  for  a  coordinate  scale 
change.  The  scale  changes  can  be  made  by  a  proportional  change  in 
various  pairs  of  variables,  where  all  other  variables  are  constant,  as 
in  the  following  algebraic  expressions. 


N* 

N  ' 

r^x  j 

N 

N*  r 

p 

Kp*  / 

when 


ap 


^  and  P  =  constant. 


X 


where 


N 

(Tp 


K,, 


=  number  of  bomblets  or  rounds  in  the  salvo 

=  delivery  error  standard  deviation  in  x-direction 

=  delivery  error  standard  deviation  in  y-direction 

=  bomblet  standard  deviation  in  x-direction 

=  bomblet  standard  deviation  in  y-direction 

=  probability  that  a  single  random  hit  on  the  target  will 
result  in  a  kill 

=  over-all  kill  probability  of  a  delivered  cluster  weapon 


^U.S.  Naval  Ordnance  Test  Station.  A  Handbook  on  the  Effective¬ 
ness  of  Cluster  Weapons  Against  Unitary  Targets,  by  Eldon  L.  Dunn. 

China  Lake,  Calif.,  NOTS,  9  December  1959.  (NAVORD  Report  7019, 
NOTS  TP  2382),  UNCLASSIFIED. 
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^F. 

F 


denotes  a  different  set  of  values  for  N,  Kp,  etc. 

ratio  indicates  the  extent  to  which  the  cluster  pattern  is 
elliptical  in  the  ground  plane. 


The  principle  obstacle  in  using  these  equations  is  the  unknown 
limits  to  which  they  will  hold. 


Several  attempts  were  made  to  explain  these  relationships  on  a 
rigorous  theoretical  basis  without  any  complete  success  (this  approach 
is  receiving  further  attention  at  NOTS). 


Eventually,  a  more  complete  expression  was  developed  v/hich 
makes  a  rather  versatile  extrapolating  tool  when  the  user  is  aware  of 

CTo 

its  limitations.  When  P  and  ■  are  held  constant, 

C^Fx 

Kp  N  L/W  Kp*  N*  L/W* 

^F^  crp  ”  CTF  *  ap  *  ' 

X  ^  y  X  ^  y 


where  L  =  target  length,  W  =  target  width  and  cr  is  expressed  in  target 
widths.  In  other  words,  when 


OTFx 


is  constant, 


I  Kp  N  L/W  \ 

\  ^  ^  / 


at  least  over  a  limited  region. 


LIMITED  THEORETICAL  ASPECTS 


Although  a  complete  theoretical  analysis  using  the  basic  equation 
is  not  available  at  the  present  time,  a  few  hints  on  the  probable  origin 
of  this  parameter  grouping  will  be  useful  in  later  discussions  of  the 
limitations  in  using  the  extrapolation  formula. 

Fortunately  it  is  possible  to  prove  that,  when  everything  else  is 
constant, 


Kp  N  =  Kp*  N* 

is  an  excellent  approximation  under  certain  conditions. 
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Fortunately  it  is  possible  to  prove  that,  when  everything  else  is 
constant, 


KpN  =  Kp«N« 

is  an  excellent  approximation  under  certain  conditions. 


If  the  expression 


|l  -  [1  -P{x.y)]Nj. 


in  the  basic  equation  (see  Appendix)  is  expanded  by  the  binominal 
series, 


(l-x)n  =  1  -  nx  + 


n(n-l)x^ 


n(n-l)(n-Z)x 


!-l)^nIxJ- 


(n-i)li! 


and  |l  -  [  1  -  KpPh  ]  Nj 


=  NKpPh  - 


N(N-l){KpPh)®  N(N-l)(N-2)(KpPh)= 


=  NK_Ph  1 


{N-l){KpPh)  (N-l)(N-2)(KpPh)‘ 


where  P^  is  the  single  round  hit  probability  defined  by 


B  u^ 


e  du 


.D  v“ 


2  j 

e  dv 


If  the  expansion  in  the  brackets,  which  is  multiplied  by  NKpPh. 
would  cancel  in  a  proportion  such  that  the  basic  equation  reduces  to 


NKp 

(2ff)® 


^  -CO  -CO  V.  L  f 
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the  approximation  that 


KpN  =  ICp*  m, 


when  everything  else  is  held  constant,  is  proven. 

That  this  expansion  in  the  brackets  does  very  nearly  cancel  in  a 
proportion  under  certain  conditions  is  shown  by  the  following  numerical 
examples. 

If  N=l,  it  is  obvious  that 

|l.  -  [1  -KpPhlNj.  =KpPh. 

If  N  is  small,  but  larger  than  1,  the  approximation  is  quite 
inaccurate. 

For  instance,  if  N  =  2,  Kp  =  0.  8,  Ph  =  0.  2  and  N*  =  4,  Kp*  =  0.4 
Ph*  =  0.2, 


{~1  -  [  1  ~  (0.8)  (0.2)  ]^y  0.2944  _  2(0.8) 
{l  -  [  1  -  (0.4)  (0.2)  ]■*}  0.  1836  4(0.4) 


.  u  J  () 


Obviously  the  ratio  of  the  bracketed  terms  is  (1.  037)  instead  of 

KpN 

1;  so  the  ratio-— =  constant  is  not  true  in  this  region. 

xvp^N^ 

However,  if 

N  -  50,  Kp  =  0.8,  Ph  =  0.2 
andN*  =  100,  Kp*  =  0.4,  Ph*  =  0.2, 


{l  -  [  1  -  (0.8)  (0.2)  ^  0.9998363  ^  50(0.8)  000075). 

{l  -  [  1  -  (0.4)  (0.2)  0.9997608  100(0.4) 

The  error  of  the  approximation  is  only  0.  0075%. 

If  N  =.  50,  Kp  =  p.  8,  Ph  =  0.  2  and  N*  =  1000,  K^,*  =  0.  04, 
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then 

-  C  1  -  (0.8)  (0.2)  ^  0.9998363  _  50(0.8)  (1.0OOI6) 

{1  -  [  1  -(0.04)(0.21  0.9996751  (1000)  (0.04) 

The  error  is  still  only  0.016%. 

If  N  =  50,  Kp  =  1,  Ph  =  1 

andN#  =  100,  Kp*  =  0.5,  Ph*  =  1, 

then 

{1  -  [  1  -  (1)(1)  _  1.00000  ^  50  (1) 

{1  -[  1  -  0.5(1)  1  -  7.8886(10“®^)  100(0.5) 

The  accuracy  of  the  approximation  is  obviously  excellent  when  the 
initial  KpPh  “•  1 

If  N  =  50,  Kp  =  0.  02,  Ph  =  0.  01 
and  N*  =  100,  Kp*  =  0.01,  Ph*  =  0.01, 

then 

{1  -  [  1  -  (0.02)(0.01)  f°}  _  0.009955  _  50(0.02)  (1.001). 

{1  -  [  1  -  (0.01)(0.01)  0.00994  100(0.01) 

Here  again,  the  approximation  error  is  small.  An  error  of  approxi¬ 
mately  0.  1%  occurs. 

The  accuracy  of  the  approximation  is  also  demonstrated  in 
Fig.  1,  2  and  3  and  will  be  discussed  further.  It  is  concluded  from 
the  preceding  examples  that  the  starting  (lower)  value  of  N  should  be. 
greater  than  50  and  the  starting  value  of  the  product  KpP^  or  of  the 
Kp 

quotient  should  not  be  extremely  small. 

Since  the  validity  of  the  parameter  grouping  (NKp)  in  the 


J 
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has  been  proven  over  the  region  of  interest,  it  is  a  little  easier  to 
speculate  as  to  the  origin  and  validity  of  the  other  parameters  in 


Kp  N  (L/W) 
ap  crp 


Kp»  N*  (L/v/)* 


^  X 


ap  * 

y 


It  is  not  difficult  to  see  the  origin  of  {crFxCrFy)  if  the  integral  is 
regarded  as  a  summation  of  an  infinite  number  of  terms--each  con- 

I 

taining,  among  other  things,  the  parameter  grouping  — 

Fx  y 

The  degree  to  which  the  remaining  terms  will  tend  to  destroy  the  pro¬ 
portion  cannot  be  shown  except  by  the  graphical  examples  that  follow 
later,  but  it  is  expected  that  for  large  values  of  OR.  and  Op,  the 
exponential  terms  in  the  infinite  expansion  will  tend  to  approach  unity, 
and  the  ratio  of  the  expansions  for  slightly  different  values  of  or..  Op 
will  be  nearly  1. 


The  inclusion  of  the  term  LW  or  L/W  when  cr  is  expressed  in 
terms  of  target  widths  is  due  mainly  to  intuitive  considerations,  It 
is  meaningless  to  speak  of  kill  probabilities  and  delivery  errors 
without  reference  to  the  target  size.  The  expression  of  delivery 
error  in  terms  of  target  width  would  necessitate  the  inclusion  of  a 
term  such  as  L/W  in  order  to  make  the  expression  dimensionally 
consistent. 


The  adjustment  of  the  target  shape,  (L/W),  can  be  expected  to  be' 
quite  limited  and  heavily  dependent  on  the  starting  values  of  CTR^^.,  crR,y. 

(or  on  F-x.>  ^Fv*  and  are  constant). 

^Fx  ^Fx 


The  correspondence  between  the  delivery  standard  deviation  and 
the  cluster  pattern  size  or  ammunition  standard  deviation  is  unique, 
as  is  the  correspondence  between  the  target  size  and  ammunition 


CTRx 

standard  deviations.  The  ratio  — — is  a  parameter  in  its  own  right, 
and  the  replacement  of  (7  r„  with  — ^  '  is  one  of  the  most  important 

^  oj?  X 


steps  in  permitting  the  handbook  graphs  to  be  superimposed  with  a 
good  fit. 


GRAPHICAL  TESTS  OF  THE  LIMITS  ON  THE 
EXTRAPOLATION  FORMULA 


The  first  test  of  proportional  relationship  between  the  various 
parameters  in  the  extrapolation  formula  is  made  with  N  and  Kp .  The 
originally  presented  extrapolation  formula  is  solved  for  Kp*. 
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were  selected  so  that 


Nl  =  10,  Kpi  =  1.  0, 

N2  =  100,  Kp2  =  0.  1, 

N*  =  100,  =  0.  1  . 

A  check  on  one  curve  will  thus  be  given. 

The  values  of  o^Fx  ~  and  L/W  =  3  were  selected  more  or  less 
arbitrarily  as  a  start. 


From  Fig.  1,  it  is  obvious  that  the  two  curves  that  were  supposed 
to  match  do  not  match  exactly.  The  curves  have  a  similar  shape,  but 
are  displaced  so  that  a  nearly  constant  error  in  P  occurs  for  the 

various  values  of-^- — . 

OtFx 


Figure  2  is  plotted  the  same  as  Fig.  1,  except  that  the  lower  value 
on  the  number  of  bomblets  is  changed  from  N  =  10  to  N  =  50,  and  Kp 
accordingly  changed  from  Kp  =  1. 0  to  Kp  =  0.  2.  For  this  plot,  there 
are  three  curves  which  can  be  used  as  a  check  on  the  extrapolation 
equation. 


From  Fig.  2,  it  can  be  seen  that  the  curves  that  are  supposed  to 
match  do  come  very  close  to  matching  over  the  entire  range  of 


O^Rx 

aFx 


(there  are  without  a  doubt  some  plotting  errors  in  the  original 


data) . 


To  be  sure  that  low  values  of  CTFx’  ^Fy  have  an  unexpected 

effect  on  the  accuracy  of  the  extrapolating  proportion.  Fig.  3  was 
plotted.  Figure  3  was  plotted  the  same  as  Fig.  2  ex.cept  that  CFx 
was  reduced  from  IW  to  0.5W.  No  degradation  in  accuracy  can  be 
detected  in  Fig.  3  as  a  result  of  this  change. 

It  is  concluded  that  as  far  as  can  be  determined  from  graphical 
tests  the  previous  numerical  evaluation  is  correct.  In  short,  the 
extrapolating  proportion 

g 

Unfortunately,  there  are  only  a  limited  number  of  handbook 
parameter  combinations  which  can  be  used  for  checking  the  extrapola¬ 
tion  formula.  Graphs  that  are  mpre  appropriate  could  be  constructed 
by  interpolating  the  handbook  data;  however,  this  would  uitroduce  an 
error  of  unknown  magnitude  which  would  be  undesirable  for  the  present 
purpose. 
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O-Fx 


FIG.  3.  A  Test  of  the  Proportional  Relationship  Between  N  and  Kp 
for  CTFx  =  o-Fx*  =  0'5W,  L/W  =  (L/W)*  =  3,  OFy/c^ Fx  =  1- 5,  When 
N  is  Large  (N  =  50)  and  is  Small  (0.  5W). 

Figure  4  also  shows  the  importance  of  the  starting  value  of  i 

that  is  chosen  when  extrapolating  with  CfFx*  region  of  interest  | 

is  to  the  right  of  the  optimum  value  of  P  for  the  various  values  of  Kp, 
the  starting  value  is  not  too  critical  and  ZW  is  apparently  adequate. 

If  the  region  of  interest  is  to  the  left  of  optimum  P,  the  starting  point 
is  very  important.  For  small  extrapolation  intervals,  a  small  value 
of  aPx  can  be  used:  but  if  a  large  extrapolation  is  made,  the  initial 
value  of  CTFx  should  be  approximately  8W.  It  is  recommended  that 
Fig.  4  be  used  as  a  guide  to  point  out  the  regions  of  limited  extrapola¬ 
tion  accuracy. 

It  is  to  be  noted  that  the  data  of  Fig.  4  for  CtFx  "  ®W  and  N  =  1600  ’ 

are  not  found  in  the  handbook,  because  these  data  are  the  result  of  . 

running  the  IBM  Program  of  the  basic  equation  for  a  check  on  the  t 

extrapolation  technique.  f 

•  S 


R  SALVO  KILL  PROBABILITY 


FIG.  4.  Reduction  of  Three  Sets  of  Kill  Probability  Curve 
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Figures  5-8  are  intended  to  show  the  limitations  in  using  the 
extrapolation  formula  to  account  for  the  change  in  the  target  shape. 

Figure  5  shows  the  handbook  graph  of  N  =  30,  L/W  =  1  (target  1) 
superimposed  on  the  hahdbook  graph  ofN=10,  L/W=3  (target  3)  for 

=  IW  in  both  cases.  The  canceling  effect  of  N  and  L/W  in  the 
extrapolation  formula  causes  Kp*  to  equal  Kp.  The  match  is  not 
especially  good  for  this  trial  value  of  ^  constx  acted. 

Figure  6  has  identical  values  of  all  parameters  except  for  (7Fx» 
which  is  increased  to  2W.  This  time  the  match  is  good  except  where 
CR 

- 2L  is  smaller  than  that  for  optimum  P. 

Figures  7  and  8  are  similar  to  Fig.  "5  and  6  except  that  the  target 
length  of  one  set  is  increased  to  5W  when  N  =  10,  and  this  is  compared 
to  a  set  which  has  a  target  length  of  IW  and  N  =  50. 

The  match  for  Fig.  7  is  very  poor,  and  for  Fig,  8  the  match  is 
fair  over  the  portion  to  the  right  of  the  optimum  value  of  P, 

It  is  concluded  that  a  good  deal  of  caution  should  be  exercised  in 
extrapolating  to  different  target  shapes.  The  value  of  ff  Fx  should  be 
as  large  as  possible,  and  each  problem  should  be  checked  separately 
for  the  error  in  making  the  extrapolation. 

The  handbook  gives  a  method  for  changing  the  target  shape  or  the 
dive  angle;  however,  it  would  have  the  same  limitations  as  have  just 
been  discussed  because  the  two  techniques  are  actually  identical. 


By  changing  the  symbol  definitions 

(L/W)*  _  °^Fx* 

L/W  CTFx  CTFy 

to  (L/W)*  =  K  (1,  3  or  5  target  widths), 

L/W  =  L  (actual  length  of  target  in  units  of  target  width), 

where  6  =  desired  dive  angle, 

or  CSC  6  for  the  handbook  dive  angle, 

X 


CTF, 


JL  - 


=  CSC  6 


and 


CfFy* 
or  T?  «  ■ 


=  1.5 
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0.8 


0.6 


0.4 


0,2 


o'Rx 

OTFx 


FIG.  8.  A  Test  of  the  Proportional  Relationship  Between  L/W  and  N 

for  0-  Fx  =  Of  Fx*  =  OfFy  =  ty  _  =  1.5,  Kp  =  Kp*. 


CONCLUSIONS 


The  relationship 

Kp  N  L/W  _  Kp*N*{L/W)« 

O'Fx^'Fy  CTFx^OTFy* 

Of  R-x 

when  P  and  — — —  c^re  constant,  can  be  used  for  extrapolating  the  data 

cr  Fx  ■ 

of  NAVORD  Report  7019  quite  accurately,  if  the  limitations  which  are 
discussed  in  the  text  of  this  report  are  well  understood. 
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Appendix 

BASIC  EQUATION 


The  basic  equation,  taken  from  the  referenced  Handbook  (NAVORD 
Report  7019)1  is  as  follows: 


where 


-bv\®  /y-b/®' 


(TF. 


D 


•'  S' 
e  dv 


and 


N 

C'F, 

b,: 

b.. 


B 

C 

D 

(yRx 

O^Ry 

X: 

Yi 

XL 

YL 


-  Xi  -X 


CTR 
=  A  + 


X 


XL 

O'Rv 


-  Yi  -y 

<^Ry 

=  c  + 


YL 


=  round  standard  deviation  in  x-direction 
=  round  standard  deviation  in  y-direction 
=  x-coordinate  of  lower  left  corner  of  the  target 
=  y-coordinate  of  lower  left  corner  of  the  target 
=  length  of  the  target  in  the  x-direction 
=  length  of  the  target  in  the  y-direction. 


=  number  of  rounds  in  the  salvo 
=  delivery  error  standard  deviation  in  x-direction 
=  delivery  error  standard  deviation  in  y-direction 
=  x-coordinate  of  the  aim  point  of  the  salvo 
=  y-coordinate  of  the  aim  point  of  the  salvo 


dx  dy 
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limitations  in  using  the  formula,  however;  and  limitations  in  using  the  formula,  however;  and 

since  tne  formula  has  an  empirical  origin,  its  since  the  formula  has  an  empirical  origin,  its 

limitations  are  discussed  at  l.ngth  with  the  aid  of  limitations  are  discussed  at  length  with  the  aid  of 

some  theory,  several  graphs,  and  some  numerical  some  theory,  several  graphs,  and  some  numerical 

example  s .  example  s . 


